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The  addition  of  artificial  viscosity  terms  to  the  Euler  equations  not  only  ensures  nu¬ 
merical  convergence,  but  also  models  the  physical  problem  correctly.  Specifically,  the 
additional  terms  guarantee  that  an  entropy  inequality  is  satisfied.  This  can  be  achieved 
by  upwind  differencing  in  the  supersonic  region  of  a  transonic  flow  so  that  the  truncation 
error  represents  an  artificial  viscosity  satisfying  the  entropy  inequality.  Most  Euler  solvers, 
however,  employ  central  differences  for  aU  spatial  derivatives,  and  therefore  expUcit  arti¬ 
ficial  viscosity  terms  are  added  to  all  four  equations  in  order  to  guarantee  convergence. 
A  check  on  the  physical  validity  of  the  additional  terms  is  to  show  that  not  only  is  the 
entropy  inequality  satisfied  for  the  modified  equations,  but  it  also  results  in  a  formula  for 
the  wave  drag  depending  on  the  specific  form  of  the  artificial  viscosity. 

Both  potential  and  Euler  solvers  in  our  transonic  aerodynamics  codes  make  use  of 
the  entropy  inequality  to  measure  the  wave  drag.  This  calculation  involves  summing  a 
positive  definite  quantity  over  the  region  of  fiow  in  which  the  shock  is  smeared,  avoiding 
the  use  of  computed  flow  variables  near  the  trailing  edge  where  there  is  uncertainty.  In 
more  standard  approaches  the  pressure  is  integrated  around  the  airfoil  or  over  a  contour 
around  the  shock,  but  the  double  integration  is  more  accurate. 

Consider  the  modified  Euler  equations 

Pt  +  +  (^)»  = 


{fm)t  -1-  (pu^  +  p)r  +  {p^v)y  =  ^ 

{pv)t  -f-  (puv)x  +  {pv^  +  p),  =  V  •  uV{pv) 
{pE)t  +  {puH)i  -f  {pvH)f  =  V  •  vV{pH)  , 


where  artificial  viscosity  has  been  added  on  the  right  of  the  form 
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These  equations  may  be  combined  to  form  a  fifth  equation  for  the  entropy 

{pS)t  +  (puS)^  +  (pvS),  =  ^  {V  •  uV(pH)  -{E-  TS)V  •  i/Vp 

-u(V  •  vV{pu)  -  uV  •  vVp)  -  u(V  •  vV{pv)  -  vV  •  vVp))  , 

whose  left-hand  side  is  in  conservative  form.  The  right-hand  side,  however,  may  be  sep¬ 
arated  into  a  positive  definite  term  and  a  divergence  term  leading  to  wave  drag  norms  of 

the  form 


{pS)t  +  {^pS)x  +  {vpS),  =  \\u  -b  IIA  A17||2  +  V  .  [uSVp  +■■■). 

An  application  of  the  divergence  theorem  then  3delds  a  representation 

Cd  =  //  (ll-VPlp +  ||AA:;|P)<ii<iy 

for  the  wave  drag  as  an  integral  of  positive  terms  over  any  region  including  the  shocks.  A 
coroUary  of  this  analysis  is  that  the  numerical  solution  of  the  Euler  equations  has  to  satisfy 
the  entropy  inequality  asserting  that  5  increases  through  shocks.  The  entropy  integral  for 
the  drag  has  the  advantage  in  three-dimensional  fiow  that  it  enables  one  to  distingmsh 
the  wave  drag  from  the  induced  drag,  which  may  be  much  larger.  On  a  practical  mesh  it 
turns  out  that  the  entropy  formula  pro-vddes  a  good  estimate  of  drag  in  the  realistic  range 
of  Mach  numbers. 

Methods  for  the  design  of  supercritical  wing  sections  have  led  to  lifting  airfoils  for 
which  two  distinct  solutions  of  the  Euler  equations  can  be  found  with  visibly  different 
shocks  at  identical  flow  conditions.  Recently  we  have  applied  complex  characteristics  to 
construct  a  symmetric  shockless  airfoil  that  has  lift  at  zero  angle  of  attack  caused  by 
shocks  on  the  upper  and  lower  surfaces  that  are  not  the  same.  Here  agziin  the  solution 
of  the  Euler  equations  is  not  unique.  The  large  size  of  the  supersonic  zones  swept  out  by 
the  characteristics  is  presumably  what  accounts  for  the  existence  of  multiple  solutions  at 
off-design  conditions. 
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